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ABSTRACT 



We report on the discovery by the Swift Gamma-Ray Burst Explorer 
of the eighth known transient accretion-powered milhsecond pulsar, SWIFT 
J1756.9— 2508, as part of routine observations with the Swift Burst Alert Tele- 
scope hard X-ray transient monitor. The pulsar was subsequently observed by 
both the X-Ray Telescope on Swift and the Rossi X-Ray Timing Explorer Pro- 
portional Counter Array. It has a spin frequency of 182 Hz (5.5 ms) and an 
orbital period of 54.7 minutes. The minimum companion mass is between 0.0067 
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and 0.0086 M©, depending on the mass of the neutron star, and the upper hmit 
on the mass is 0.030 Mq (95% confidence level). Such a low mass is inconsistent 
with brown dwarf models, and comparison with white dwarf models suggests that 
the companion is a He-dominated donor whose thermal cooling has been at least 
modestly slowed by irradiation from the accretion flux. No X-ray bursts, dips, 
eclipses or quasi-periodic oscillations were detected. The current outburst lasted 
~ 13 days and no earlier outbursts were found in archival data. 

Subject headings: binaries: close — pulsars: general — stars: neutron — white 
dwarfs — X-rays:binaries 



Introduction 



It is now well established that a subset of neutron stars in low-mass X-ray binary systems 
are millisecond pulsars powered by accretion from either a brown dwarf or degenerate com- 
panion in a very close bi nary conflguration . Such systems are believed to be progenitors of 
millisecond radio pulsars (lAlpar et al.lll982l ). Before the discovery of SWIFT J1756.9-- 2508, 
there were seven known accretion-powered millisecond pulsars (see lKaaret et al.ll2006l . for a 
full listing). The sample is now becoming large enough to start to study classes of these ob- 
jects which can in turn lead to a better understanding of the properties of neutron star X-ray 
binaries and their relationship to millisecond radio pulsars. The known accretion-powered 
millisecond pulsars all appear to be in very close orbits (with orbital pe riods ranging between 



41.1 and 256.5 min). For all but one such system, XTE J1814— 338 (IBhattacharyya et al. 



20051 ). orbital constraints rule out a hydrogen-burning main sequence companion, and for 
four systems, even a brown dwarf companion is not consistent with the orbital parame- 
ters. The application of recent models of w hite dwarf donors in ultracompact binaries (e.g. 
Deloye fc Bildstenll2003l : iDeloye et al.l 120071 ). shows that, within certain temperature ranges, 
white dwarfs are a viable model for the companion in the most compact of the X-ray mil- 
lisecond pulsar systems. 

In this Letter, we report the discovery of SWIFT J1756.9— 2508 (hereafter referred to 
as J1756), with a frequency of 182.07 Hz. It was discovered by the Swift Gamma-Ray Burst 
Explorer Burst Alert Telescope (BAT) as part of the hard X-ray transient monitor program. 
In § [2] we discuss the observations made with Swift and the Rossi X-ray Timing Explorer 
(RXTE). In § 13] and § H] we discuss the timing and spectral analysis and in § O we derive 
limits on the companion mass and compare this source to the other known accreting X-ray 
millisecond pulsars. All errors are at the 90%-confldence level unless otherwise stated. 
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2. Observations and data reduction 
2.1. Swift 

SWIFT J1756.9 — 2508 was discovered using the SwiftBAT hard X-ray transient monitor 



(iKrimm et al.ll2006l ). The BAT is a large field of view instrument that continually monitors 
the sky to provide the gamma-ray burst trigger for Swift. On average more than 70% of 
the sky is observed on a daily basis. Sky images are processed in near real-time to detect 
astrophysical sources in the 15-50 keV energy banc0. 

J1756 was detected in the daily full-sky mosaics that are produced in the transient 
monitor. It first rose to detectability (Fig. [1]) on 2007 June 7 (MJD 54258) with an aver- 
age rate in Swift BAT of 0.0083 ± 0.0015 ct c m~^, corresponding to approximately 37 



mCrab. The discovery was announced on June 13 (IKrimm et al.l 120071 ). The BAT data for 
all pointings when the source was fully coded (visible to all BAT detectors) were processed 
using the Swift analysis tools, which included cleaning out the diffuse background and flux 
from other bright sources in the fleld of view. 

The discovery of J1756 triggered a target of opportunity (ToO) observation in the Swift 
X-ray Telescope (XRT) beginning at 15:25:32 UT, 2007 June 13. It was intermittently ob- 
served over 13 days for a total exposure of 60 ks. The XRT data were processed with standard 
procedures (xrtpipeline vO.10.6 with in FTOOLS in the H easof t package v. 6. 1.2), flltering. 



screening, and grade selection criteria (IBurrows et al.ll2005l ). The photon counting (PC) data 



collected during the flrst three segments were corrected for pile-up and for the presence of 
single-reflection rings due to the nearby bright source GX 5-10 For the windowed timing 
data we extracted source events in a square region with a side of 20 pixels. Ancillary response 
flies were generated with xrtmkarf and account for different extraction regions, vignetting 
and Point-Spread Function (PSF) corrections. The light curve (Fig. [T]) was corrected for 
PSF losses and for vignetting. The best flt XRT position is shown in Table [TJ 



2.2. RXTE 

The discovery of J1756 also triggered a ToO observation in the Rossi X-ray Timing 
Explorer (RXTE), starting at 18:42:00 UT, 2007 June 13. Because GX 5-1 is within the 
^ 1° fleld of view of the RXTE Proportional Counter Array (PCA) collimator, the PGA 



'See http://swift.gsfc.nasa.gov/docs/swift/results/transients/ for details on the transient monitor. 
'XRT observations after MJD 54278 are not usable due to severe contamination from GX 5—1. 
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pointing direction was offset by about 18', primarily eastward of the true pulsar position. The 
PCA light curve (Fig. [1]) has been corrected to effective on-axis values using the collimator 
response function, and spectral response functions account for the offset as well. The offset 
pointing direction is within 3' of the Galactic ridge, and so contamination from diffuse 
emission is present. The amount of contamination was found by examining other nearby 
PCA observations, when J1756 was quiescent, and fitting with a simple phenomenological 
model (and assuming it is constant). The contaminating flux was about ~3 mCrab. The 
total exposure was ^104 ks in two observing programs. 

Archived monitor data from RXTE and Swift were searched for earlier outbursts of 



J 175 6. The PCA performs regular monitoring observations of the Galactic bulge (ISwank fc Markwardt 



20011 ) and a search of archived PCA bulge scan data shows that the source has not been 
detected in the 2-10 keV band over the past 8.4 years, with an upper limit of about 5 mCrab 
(2-lOkeV)0 This limit is likely to be systematics driven, based on the nearby bright source 
GX 5—1. Archived daily BAT monitor data (15-50 keV) were searched back 2.5 years and 
the i?Xri? All-Sky Monitor light curves (2-10 keV) were searched back 11.4 years with no 
evidence of the source in either instrument with Sex upper limits of ~ 15 mCrab (BAT) and 
f» 50 mCrab (ASM). Due to observing constraints (mostly for sun avoidance), the overall 
good coverage of the three detectors is ~ 87%. With these gaps and the fact that even the 
2007 outburst was not seen strongly in the ASM, it is possible that an earlier outburst was 
missed. Thus we set a tentative recurrence time scale of > 10 yr. 



2.3. Other Observations 



The discovery announcement led to several multiwavelength observations of J1756. A 
possible fading near-IR cou nterpart was annou nced with Kg = 19.7 on MJD 54270 and 
Ks = 21.0 on MJD 54282 jBurderi et al.ll2007h . However an observation with the West- 
erbork Synthesis Radio Telescope during the X-ray outburst (1.4 GHz; MJD 54266.103; 
Hessels fc Stapperd 120071 ) showed no sign of coherent radio pulsations. Furthermore, J1756 
was not detected in observations after the outburst in either 8.7 GHz radio (MJD 54277.371, 



Possenti et al.ll2007l ) or 0. 3-8.0 keV X ray (C handra upper limit of 1.4 x 10 erg cm ^ s ^, 



starting MJD 54287.059. IPapitto et alll2007h observations. 



^''The 2007 outburst of J1756 was not seen in the PCA bulge scan because RXTE spacecraft slewing 
constraints prevented observations in early June. The most recent observation before the outburst was 2007 
June 6. 
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3. Timing Analysis 



Fig. [Ushows that the outburst of J1756 had a very rapid rise (< Id). After this the BAT 
flux was nearly constant for the next several days. For the first several days of XRT and PCA 
observations (up to MJD ^ 54268.6), the flux in both instruments showed a slow decline, 
which can be fit to a nearly exponential decay, e"*/"^, where r = 7.3 ± 0.4 d. This matches 
fairly closely the decay in the BAT light curve. After this date the flux in all three instruments 
began a steeper decline (r = 0.6 ± 0.3 d) and was undetectable even in the PCA after MJD 
54271, a mere 13 days after the start of the outburst. This light c urve is similar in form to 
those of three other X-ray millis econd pulsars: SAX J 1808.4— 3658 (IBildsten fc Chakrabarty 
200lh and ICR J00291 +5934 JCallowav et all l2005h . and in particular XTE J1751-30 5 
( iMarkwardt et al.ll2002l ). which showed very similar values of r ( IGierlihski fc Putanenll2005l ). 
The three other objects have had more than one outburst, on recurrence time scales of ~ 2— 3 
yr, and all outbursts of XTE J1751-305 and ICR J00291+5934 and most o utbursts of SAX 
J180 8.4— 3658 h ave been short (< 15 d ays), although notably the 2000 (IWijnands et al. 



20011) and 2002 (IMarkwardt et al.ll2005r ) outbursts of SAX J1808.4-3658 were longer and 
more erratic. Given the non-detection of earlier outbursts, J1756 is alone among the fast- 
decaying millisecond pulsar systems in not showing recurrence. 

We searched for high frequency variability by constructing an FFT-based power spec- 
trum from the PCA da ta. The spectrum showed a significant, narrow excess near 182 Hz 
( IMarkwardt et al.l 120071 ). which was also found in following observations. The feature was 
several mHz wide, which suggested an X-ray pulsar with orbital Doppler broadening. Exami- 
nation of power spectra on short (f»300 s) intervals confirmed a sinusoidal orbital modulation, 
making SWIFT J1756.9— 2508 the eighth known accreting millisecond X-ray pulsar, and the 
one with the longest rotation period. 



m 



We performed precisio n timing analysis on the PCA data using the techniques described 



Markwardt et al. 



(120021 ). X-ray event data were recorded in the E_125us_64M_0_ls mode, 
and event times were corrected to the solar system barycenter using the position determined 
from XRT imaging (see Tabled]). The Z"^ statistic (or Rayleigh statistic) was computed, and 
the orbital/timing parameters were adjusted to produce the largest power of the pulse 
fundamental. This p rocedure is equivalent to radio or X-ray pulse time-of-arrival fitting 
( Buccheri et al.lll983l ). but with no binning. After adjustment, we found no long term or 
orbital systematic trends. The binary orbit model was the ELLl from TEMPO software 
version 11.0100 which is designed for nearly circular orbits. The orbital parameters are 
shown in Table [U and confirm that this is a very compact and circular orbit. The pulsation 



^'"^ http://www.atnf. csiro.au/research/pulsar/tcmpo/ 
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frequencies from MJD 54264.78 to 54267.07, folded on a trial orbital period clearly establish 
the orbital frequency modulation. The pulsed semi-amplitude is a nearly constant ~6% 
throughout the observations. No pulsed signal was seen in the 15-50 keV band from 250 s of 
BAT event starting starting at MJD 54266.065, with a 2a limit of < 13 %. With the present 
data we can place only weak upper limits on the pulse frequency derivative (/), since the 
time basehne is quite short. 

The data from both the PCA and XRT (which cover 100% of the orbit) were searched 
for modulation at the orbital period. The PCA light curve (2-10 keV) shows no modulation 
at the orbital period, after the long term trend is subtracted. The 2a upper limit to orbital 
modulations is < 1.2 % from the PCA data. There is also no significant modulation seen in 
the XRT data. Furthermore there are no signs of eclipses, dips or bursts. We searched for 
rapid variability by constructing power spectra of 64 s segments of event data, with a Nyquist 
frequency of 2048 Hz, and averaged contiguous power spectra. While some variability is 
present, especially below ~200 Hz, a detailed timing study is beyond the scope of this 
paper. We do not find any features above 200 Hz which are obvious kiloHertz quasi-periodic 
oscillations (QPOs). For the summed outburst, the 95% upper limits on the fractional r.m.s. 
are about 6% (9%) for 16 Hz (128 Hz) wide QPOs, and a factor of ~2 larger for a single 
observation. 



4. Spectral Analysis 

We fit the Swift/XRT spectrum from each segment in the 0.5- 9 keV energy range with 



different models, which included a single absorbed power law (PL; IWilms et al.ll2000l ). and 
an absorbed black body (BB). A simple comparison of null hypothesis probabilities clearly 
favours the PL with respect to the BB model (20% vs. 4%, respectively), with the black 
body fit producing systematic residuals both at low and high energies. The derived ab- 
sorbing column density i s always in excess of the Galactic value (A^^ = 1.47 x 10^^ cm~^. 



Dickey fc LockmanI Il990l ) and there is no indication of a softening of the spectrum with 
time. Representative spectral fits are shown in Table [2l The source was detected in the 
PCA up to energies ~ 30 keV, and the spectrum is consistent with an absorbed power law 
with absorption fixed at the best-fit value of XRT segment 003. The BAT data from 14 to 
195 keV were fit by a simple power law model. The BAT spectrum is softer than the spectra 
found from fits to the PCA alone or XRT alone, suggestive of a cut-off to the power law. 
However joint fits to the BAT, XRT and PCA spectra are not improved with a cut-off power 
law model. 



Given the uncertainty in the spectrum at high energies, we estimate the integrated 
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outburst fluence to be ~ (4.5 ± 0.8) x 10~^ erg cm~^ (1-10000 keV). Using the angular 
proximity of the source t o the Galactic center to estimate the distance, we can, following 
Bildsten fc Chakrabartyl (1200 ll ). estimate the time-averaged mass accretion rate, Mx = 
(9.3 ±1.6) X 1O~^^M0 jY'^dlm^^T^^ , where dg is the distance is units of 8 kpc, mi. 4 = 
Mx/(1.4 Mq), and we assume a neutron star radius of 10 km. The parameter Tio is the 
recurrence time in units 10 yr, and only Tiq > 1 is allowed by the observations. 



5. Discussion 



In combination, the measured mass function and binary orbital period Pb constrain the 
donor's mass-radius {Mq-Rc) relation, assuming the donor fills its Roche lobe. For SWIFT 
J1756.9— 2508, Pb = 54.7 minutes, a value near the median of the known X-ray millisecond 
pulsars. If we assume a particular value for Afj;, we can calculate possible values of Mc from 
the mass function, = {McSiniY /{M^ + M^)"^, where i is the binary inclination to 
our line of sight and and M^. are the masses of t he pulsar arid com panion, respectively. 
The lack of eclipses allows us to set a limit i < 85° f lPaczyhskil Il97ll ) . For this inclination 
and Mx = 1.4M0, we have Mc^mm = O.OO67M0. Assuming a uniform distribution of 
possible values of cosz we can set a 95% confidence level (CL.) upper limit on the mass of 
Mc < O.O22M0 (for Mx = I.4M0). The range of possible values of M^ up to this limit 
and corresponding companion radius Rc are shown in Fig. [2J If we consider a larger neutron 
star of Mx = 2.2Mq, the limits (95% C.L.) are M^ < O.O3OM0 and Rc < 0M9Rq. This 



radius is less than the minimum radius for low-mass 



age (jChabrier et al.ll2000l: iBildsten &: Chakrabartv 



2001 



lydrogen rich brown dwarfs of any 
, Figure 2), as ex pected for a binary 



period of < 80 min ( Paczyhski fc SienkiwwiczT 1981 : Nelson et al. 1986 ). Thus we restrict 



our attention to white dwarf (WD) donor models. 

The dotted and dashed curves in Fig. [2] show evolutionary sequences of ultracompact 
binaries with He donor models originally developed for a study of the AM CVn (double WD) 
class of ultracompact binaries (IDeloye et al.l 120071 ). Along these curves, regions of negative 
slope correspond to adiabatic donor evolution. In regions of positive slope the donor is 
cooling and contracting towards the fully-degenerate He WD Mc-Rc relation (see Deloye 
et al., 2007), essentially traced by the lowest dotted line in Fig. [2l The donor's cooling 
is affected by the (very uncertain) efficiency with which irradiating flux generated by the 
accretion — which always dominates the donor's own flux in these ultracompact binaries — is 
thermalized in the donor's upper atmosphere. The impact of irradiation is to extend the 
donor's adiabatic evolution phase and to slow its rate of cooling afterwards. 



Since calculations similar to the 



Deloye et al, 



(I2OO7I ) He- donor study do not yet ex- 
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ist for C/0 donor s , we compare the J1756 constraints to a sample of pure C models of 



Deloye &: BildstenI (120031 . solid curves in Fig. [2]); realistic C/0 WDs will lie at somewhat 
smaller Rc- These models assume fuUy-convective interior profiles and do not explicitly 
treat thermal transport. However, the physics that produces the transition from adiabatic 
evolution to cooling in He donors does not depend on composition and C/0 donors will also 
experience a cooling phase. This means a substantial fraction of the phase space for "hot" 
C/0 donors shown in Fig[2]will likely not be evolutionarily accessible (e.g., the nearly vertical 
regions of the log(Tc/K) > 6.4 "isotherms"), making a C/0 donor in J1756 rather unlikely. 
Thus, J1756 likely harbors a He-dominated donor whose thermal cooling has been at least 
modestly slowed. Irradiation provides a natural mechanism that is able to affect the neces- 
sary degree of slowing, as discussed above. This interpretation places SWIFT J1756.9— 2508 
in the company of at lea st two other X-ray mi llisecond pulsars with ve ry low mass compan - 



ions: XTE J0929-314, flGallowav et all 120021 ) and XTE J1807-294 flFalanga et al.ll2005[ ): 



with its longer orbital period, J1756 provides a unique probe amongst the known ultracom- 
pact LMXBs and appears to indicate that irradiation plays a least a modest role in the 
evolution of the donors in these systems. 
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grant number TM7-800. We acknowledge Jean Swank and the RXTE Science Operations 
Facility and the Swift mission operations team for their help in rapidly scheduling these 
target-of-opportunity observations. 



REFERENCES 

Alpar, M. A., Cheng, A. F., Ruderman, M. A., & Shahm, J. 1982, Nature, 300, 728 

Bhattacharyya, S., Strohmayer, T. E., Miller, M. C, & Markwardt, C. B. 2005, ApJ, 619, 
483 

Bildsten, L. & Chakrabarty, D. 2001, ApJ, 557,292 

Buccheri, R., et al. 1983, A&A, 128, 245 

Burderi, L. et al. 2007, Astr. Tel., 1132 

Burrows, D. N. et al. 2005a, Space Science Review, 120, 165 



-9- 

Chabrier, G., Baraffe, I., AUard, F., & Hauschildt, P. 2000, ApJ542, 464 
Deloye, C. J. & Bildsten, L. 2003, ApJ, 598, 1217 

Deloye, C. J., Taam, R.E., Winisdoerffer, C, & Chabrier, G. 2007, MNRAS, accepted (astro- 
ph/0708.0220) 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

Falanga, M. et al. 2005, A & A, 436, 647 

Galloway, D. K., Chakrabarty, D., Morgan, E. H., & Remillard, R. A. 2002, ApJ, 576, L137 

Galloway, D. K., Markwardt, C. B., Morgan, E. H., Chakrabarty, D. & Strohmayer, T. E. 
2005, ApJ, 622, L45 

Gierlinski, M., & Putanen, J. 2005, MNRAS, 359, 1261 

Hessels, J. W. T. & Stappers, B. W. 2007, Astr. Tel., 1129 

Kaaret, P., Morgan, E. H., Vanderspek, R., & Tomsick, J. A. 2006, ApJ, 638, 963 
Krimm, H. A. et al. 2006, Astr. Tel., 904 
Krimm, H. A. et al. 2007, Astr. Tel., 1105 

Markwardt, C. B., Swank, J. H., Strohmayer, T. E., in't Zand, J., & Marshall, F. E. 2002, 
ApJ, 575, L21 

Markwardt, C., Swank, J., Wijnands, R., & in't Zand, J. 2005, Astr. Tel., 505 

Markwardt, C., Swank, J. H., & Krimm, H. A. 2007, Astr. Tel., 1108 

Nelson, L. A., Rappaport, S. A., & Joss, P. C. 1986, ApJ, 304, 231 

Paczyhski, B. 1971, ARA&A, 9, 183 

Paczyhski, B. & Sienkiwwicz, R. 1981, ApJ, 248, L27 

Papitto, A. et al. 2007, Astr. Tel., 1133 

Possenti, A. et al. 2007, Astr. Tel., 1128 

Swank, J., & Markwardt, C. 2001 in ASP Conf. Ser. 251, New Century of X-Ray Astronomy, 
ed. H. Inoue h H. Kunieda (San Francisco, ASP), 94 



-10- 



Wijnands, R., Mendez, M., Markwardt, C, van der Klis, M., Chakrabarty, D., & Morgan, 
E. 2001, ApJ, 560, 892 

Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 11 - 



Tabic 1. Timing Parameters of SWIFT J 1756. 9 -2508 



Parameter 



Value 



Right ascension, a (J2000.0)'' 

Declination, S (J2000.0)^ 

Barycentric pulse frequency, /o (Hz) 

Pulsar frequency derivative, |/| (Hz s"'^) 

Projected semimajor axis, sini (It-ms) 

Binary orbital period. Pi, (s) 

Time of ascending node, Tasc'' 

Orbital eccentricity, e " 

Pulsar mass function, (IQ'"^ Mq) <= 

Minimum companion mass, Mc (10^"* Mq) 

Maximum Power, Z5,„^ 



25°06'27('8 
182.065804253(72)^ 

< 1 X 10-12C 

5.942(27) 
3282.104(83) 
54265.28707(6) 
< 0.026 
1.56(3) 
6.7-9.2 f 
2061 



^The estimated position uncertainty is 3'.' 5. Position is from the 
XRT and held fixed as a timing parameter. 

'^Uncertainties are 1<t in the last quoted digits. 

^^95% upper limit. 

''Modified Julian days, referred to barycentric dynamical time. 
"^Derived parameter. 

^For neutron star masses of 1.4 — 2.2M0. 



Table 2. Sample spectral fits for SWIFT J1756.9-2508 





M.li) 


Alndc-i" 


An 


r >■:/;'' 


Arid/ 


XRT 


54266.06 


PL 


5.44 ± 0.41 


2.08 ±0.14 


1.081/200 


XRT 


54266.06 


BB 


2.82 ± 0.23 


[1.30 ± 0.05] 


1.185/200 


PCA 


54264.79 


PL 


5.44"= 


2.00 ±0.14 


0.878/47 


PCA 


54267.46 


PL 


5.44"= 


1.95 ±0.08 


0.514/48 


BAT 


54258 


PL 




2.3 ±0.7 


0.20/6 


BAT 


54266 


PL 




2.8 ±0.5 


0.38/6 



''PL=power-law; BB=absorbed black body 

''Numbers in column 5 are P for PL fits and kT (in square brackets) 
for the BB fit. 

''Parameter held fixed. 
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54255 54260 54265 54270 54275 54280 54285 

MJD 



Fig. 1. — The combined light curve of SWIFT J1756.9— 2508 from Swift BA T (diamonds). 
Swift XRT {filled circles), RXTE PCA {open circles), and Chandra (triangle: iPapitto et al 



20071 ). The BAT points show flux in the 15-50 keV band, Chandra the 0.3-8.0 keV band, 



and the other points the 2-10 keV band. The PCA data is contaminated by emission from 
the Galactic ridge. To correct for this, a constant baseline of 1.24 x 10~^° erg cm~^ s~^ 
was empirically determined and subtracted from the PCA data. Data points with < 2a 
significance are shown as upper limits. For comparison with the archival upper limits quoted 
in §2.2, for the PCA (BAT) 1 mCrab ^ 2.4 (5.2) x lO^^^ erg cm'^ s'K 
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0.005 0.010 0.015 0.020 

Companion Moss (M©) 



Fig. 2. — Companion radius (Re) vs. mass (Mc) plane, showing the Roche lobe constraints 
for SWIFT J1756.9— 2508 and, for comparison, three other X-ray millisecond pulsars with 
low-mass companions. The source curves trace possible values of the inclination angle, i. The 
solid black dots indicate i = 60° and the open black dots indicate i = 30°. The other curves 
are various white dwarf models. The d otted and dashed lin es show evolutionary sequences 
with He donor models calculated as in iDeloye et al.l (120071 ) (see text for full explanation). 
The dotted curves show cases where no accretion flux is reprocessed. The dashed curves show 
the case where 10% of accretion flux, as seen by the donor, is thermalized in the donor's 
atmosphere. The numbers next to these cu rves indicate log -^n , where higher values of the 
central degeneracy parameter ipc (defined in lDeloye et al.l 120071 ) indicate a more degenerate 
and hence compact donor. Within each set of models, evolution proceeds to lower Mq and 
less degenerate do nors evolve at larger Rg- The solid curves are for a pure C composition 
from the models of iDeloye fc BildstenI (120031 ). The number next to each solid curve indicates 
the log of the temperature in Kelvin used in the model. 



